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A theory of electronic Raman scattering in nearly antiferromagnetic Fermi liquids is constructed 
using the phenomenological electron-electron interaction introduced by Millis, Monien, and Pines. 
The role of "hot spots" and their resulting signatures in the channel dependent Raman spectra is 
highlighted, and different scaling regimes are addressed. The theory is compared to Raman spectra 
taken in the normal state of overdoped Bi2Sr2CaCu2 0g+6, and it is shown that many features of 
the symmetry dependent spectra can be explained by the theory. 



It is "widely believed that strong antiferromagnetic cor- 
relations are an important ingredient needed to describe 
the unusual properties found in the normal state of the 
cuprate superconductors. A phenomenologically based 
theory proposed by Millis, Monien, and Pines Q| sup- 
poses that the electronic system interacts with an inde- 
pendent model spin fluctuation spectrum. This nearly 
antiferromagnetic Fermi liquid (NAFL) model has been 
applied to study both spin and charge correlation prop- 
erties Of the latter, recently Pines and Stojkovic 
have investigated the two-particle correlation function in 
a perturbation approach yielding the optical conductiv- 
ity response ||. In this paper we extend their results to 
investigate the two-particle Raman response. 

The salient point of the calculation starts with assum- 
ing a form for the quasiparticle-quasiparticle interaction: 
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which at lowest order in the coupling constant g 2 yields 
for the self energy 
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with Go the bare Green's function. Here u) s f and £ are 
the phenomenological, temperature dependent spin fluc- 
tuation energy scale and the correlation length, respec- 
tively, which can be determined via fits to the magnetic 
response data B . These functional parameters obey cer- 
tain relations depending on different temperature and 
doping regimes. In the z — 1 or pseudo-scaling regime, 
the spin correlations are strong enough to lead to changes 
from the classical mean field theory z = 2 regime. As 
has been pointed out in Q, the z = 1 regime is most 
appropriate to optimally doped cuprates at intermedi- 
ate temperatures while z — 2 applies to optimally doped 
materials at higher temperatures and/or more overdoped 
materials. For the z = 1 scaling regime, uj s f£, =constant 
while for z = 2, 0J s f£, 2 = constant. Once the self energy 
is determined, in the absence of vertex corrections to the 
Raman vertex 7(k), the Raman response is simply given 
by 



X ; 7 (q = 0,tt) = 2£ 7 2 (k) f ^[/(x)-/( x + n)] 
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with G(k, x) — [x — e(k) — S(k, x)]^ 1 and band structure 
e(k) = — 2t[cos(k x a) + cos(k y a)] + At 1 cos(k x a) cos(k y a) — 
/i. As in Ref. S, we use the following approximations: 

(i) we neglect the real part of the self energy given in (2), 

(ii) we only use the Green's function evaluated at low- 
est order in the coupling constant g 2 , (iii) we neglect all 
vertex corrections for the Raman vertex, (iv) momentum 
sums are replaced by ^ k = J dQ k / | v k | J cfe(k). While 
(iv) does not crucially affect the results, approximations 
(i-iii) while simplifying the calculations considerably miss 
important quasiparticle renormalizations at larger values 
of the coupling. Therefore we expect that these calcula- 
tions would be most appropriate to describe the spectra 
taken on overdoped cuprate superconductors. We believe 
that these approximations must be lifted in order to de- 
scribe optimally and underdoped cuprate systems. This 
is analyzed in detail in a forthcoming publication Q . 

The parameters we have used for both scaling regimes 
are g = leV, a = 2.6 states/eV, t = 250meV, t'/t = 0.45, 
filling < n >= 0.8, 7Bi 9 (k) = bi[cos(k x a) — cos(A: y a)], 
7_B 2 „(k) = 62 sin(fca;a) sin(fcj,a). In addition, for the 
z = 1 scaling regime we have used uj s f^/a = 50meV 
and a^ 1 = 0.1 + 4.64T/2i, while for z = 2, uj s f£, 2 /a 2 = 
60meV and w sf /2t = 0.0237 + 0.55T/2i. These parame- 
ters are similar to those used in Ref. || to describe the 
Hall conductivity data in YBa2Cu307. For very small 
incoming laser frequencies, the Raman vertices can be 
written in terms of the curvature of the electron bands. 
However, it has been pointed out || that this is of ques- 
tionable use for the cuprates and therefore we have used 
the first order terms of a Brillouin zone expansion for the 
vertices. Therefore the magnitude for the scattering is ar- 
bitrary and determined by the dimensionlcss coefficients 
61,62, which are set by fitting to the data. This has no 
effect on the frequency dependent lineshapes however. 

Our results are summarized in Figure 1 , where we have 
plotted the B\ g and £?2 S Raman response for both scaling 
regimes. The A\ g response is slightly more complicated 
due to screening effects and is discussed in detail in Ref. 
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FIG. 1. Electronic Raman response for the B\ g 

and £>2 fl channels evaluated at different temperatures 

(T/2t = 0.01,0.02, ...,0.06) for z = 1 and z = 2 scaling 

regimes of NAFL. Here we have set a = bi = 62 = 1, as 
defined in the text. 



The spectra for both scaling regimes share several 
features. First, one immediately sees that the flat con- 
tinuum at high frequencies which is present in the Raman 
data from all cuprate superconductors as well as several 
A- 15 compounds is reproduced by the theory. This is 
a consequence of a scattering rate which is effec- 

tively linearly dependent on w at frequency scales larger 
than cjgf. Moreover, the Raman response is different for 
the different scattering geometries. This is due to the 
strong anisotropy of the scattering rate for different di- 
rections in momentum space. In the NAFL model, the 
quasiparticle scattering is strongest near "hot spots" , i.e. 
regions of the Fermi surface which can be connected by 
the antiferromagnetic wave- vector Q = (±7r, ±7r). Since 
these regions are close to the zone axes for the given 
band structure, this means that the B\ g geometry most 
effectively probes these hot spots while the B^g geome- 
try probes along the zone diagonals and therefore sees 
"colder" quasiparticles. Thus the "peak" of the spectra 
is located at larger frequencies for the B\ g channel com- 
pared to B2 g reflecting the larger energy scale. Lastly, 
the coupling constant g controls the position and charac- 
ter of the peak in each scattering geometry: for smaller 
couplings, the peak is more pronounced and is shifted to 
lower frequencies, while for larger couplings the peak is 
smeared out and pushed to larger frequencies. 

We next discuss the results for the different scaling 
regimes. For these parameter choices the differences be- 
tween the two regimes are mostly quantitative. In general 
the z = 1 results are more strongly temperature depen- 
dent than the results for z — 2. This is most clearly seen 
in the high frequency portion of each spectra and is due 
to the weaker temperature dependence of £, u> s f for z = 2 
scaling than for z = 1 scaling. In addition the z = 1 re- 
sults show stronger differences between the results for the 
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spectra for both z = 1 and z — 2 as well as the B\ g 
spectra for z = 2 show a monotonically decreasing inten- 
sity as the temperature is increased as spectral weight is 
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FIG. 2. Fit of the (z = 2) NAFL theory to the B lg spectra 
taken in g on overdoped Bi2Sr2CaCu2 0s+6 (T c = 55 A') at 
60 A, 90 A, 150 A, and 200 A. Parameters used are given in the 
text. 



being transfered to high frequencies. However, the B\ g 
spectra for z — 1 show that the low frequency portion 
decreases with the temperature but the high frequency 
portion increases with temperature. Therefore the spec- 
tral weight is being transfered out to only slightly higher 
frequencies. This may be due to the more pronounced 
"hot spots" in the z = 1 regime compared to the more 
smeared "hot spots" for z = 2 || . The B\ g channel thus 
is more sensitive to the growth of the correlation length 
as temperature is lowered than B2 g - 

As a consequence, the Raman spectra are quite sensi- 
tive to the details of the parameter choices and therefore 
the shape of the Raman spectra can be a very useful 
tool to understand and probe the anisotropic quasiparti- 
cle scattering rates in much the same way as it has been 
used to probe the anisotropy of the energy gap A(k) in 
the superconducting state 0. 

We now turn our attention to the comparison of the 
theory to the experimental data on the cuprates. As 
we remarked, the approximations used in this paper are 
most appropriate for systems with weak spin fluctua- 
tion scattering, and therefore we expect that our results 
would best represent the data from appreciably over- 
doped cuprate superconductors. Thus we consider the 
results from z = 2 scaling regime in closer detail. 

As discussed in Ref. 0, it is believed that the quasi- 
particles in the overdoped cuprates feel an effectively 
stronger impurity interaction coming from the increased 
effective dimensionality of the system which allows the 
quasiparticles to interact more strongly with defects 
residing out of the CuC>2 plane. Therefore in addi- 
tion we will consider an isotropic impurity interaction 
H im P = Ek,k'E l , CT C / e l ( k - k ') R 'cj c crCk%(T , where R, de- 
notes the position of the impurity labeled by i and U 
is the impurity potential. After averaging over the po- 
sition of the impurities, this adds a momentum inde- 
pendent term to the imaginary part of the self energy 



r imp = Trn % N F I U 



2 , where rii is the impurity concen- 
tration and Np is the density of states per spin at the 
Fermi level. 
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FIG. 3. Fit of the (z = 2) NAFL theory to the B 2g spectra 
taken in g on overdoped Bi2Sr2CaCu2 0s+<5 (T c = 55K) at 
60K, 90K, 150K, and 200K. Parameters used are given in the 
text. 



The fit of the theory to the temperature depen- 
dent spectra for each channel obtained in overdoped 
Bi 2 Sr 2 CaCu20 8+ 5 (Bi 2212, T c = 55K) by Hackl et 
al. in Ref. || is shown in Figs. 2 and 3 for the B\ g 
and B2 g channels, respectively. We note that uj s f and £ 
have not been determined via fits to magnetic response 
data. These fits are thus obtained by choosing: uj s f and 
£ at a fixed temperature, the magnitude of the impu- 
rity scattering rate, a magnitude of the coupling con- 
stant ag 2 , and &i, 62 which control the magnitude of the 
Raman scattering in each channel. These parameters 
were chosen to fit the lowest temperature data shown 
in the Figure. Then, only the temperature dependent 
parts of uj s f and £ were then modified to fit the data at 
other temperatures as shown by the rest of the curves 
in Figs. 2 and 3. The parameters we have used are 
62/61 = 0.417, T imp = 40cm- 1 ,w s/ £ 2 /a 2 = 13meV, for 
all temperatures, and uj s j = 160K + 0.06T[K]. 

Further information can be obtained concerning the 
slope of the spectra at vanishing frequencies. In Ref. ||] 
it was shown that the inverse of this slope has qualita- 
tively different behaviors for different doping regimes of 
various cuprate materials. Within the level of our ap- 
proximations, the low frequency Raman spectra is given 
by 



0) =Q 
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[2cosh(a;/2)] 2 W(k,xT) 



In Fig. 4 we plot f2/%" (fi — > 0) obtained from the 
parameters used in Figs. 2 and 3 and compare the results 
to the data taken in [g . The fit to the data is quite good 
when the error bars of the data are taken into account. A 
8% uncertainty has been included in the data || . We re- 
mark that similar fits can be made with slightly different 
parameter choices than the ones considered here. 

In closing, we remark that while good agreement be- 
tween the theory and the data has been shown for over- 
doped Bi 2212, several features remain to be explained 
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FIG. 4. Fit of the inverse slope of the Raman response for 
the two symmetry channels. 



in optimally and underdoped systems. For instance, 
an almost completely temperature independent high fre- 
quency background is observed in these systems, and the 
difference of the inverse slope of the Raman response 
between B\ g and B 2g channels grows remarkably upon 
lesser dopings ||. Moreover, features which have been 
associated with a pseudo-gap have been observed in the 
underdoped systems |l(J . These features are discussed in 
greater detail in 
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